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摘要 : 快速 射电 爆 (Fast Radio Burst 简称 FRB ) 于 2007 年 首次 在 脉冲 星 搜索 数据 中 被 发 现 ， 是 一 种 在 
射电 波段 强烈 爆发 的 天 体 物 理 现象 ， 其 特征 为 持续 时 间 一 般 仅 为 若干 毫秒 ， 峰 值 流量 密度 可 达到 央 斯 基 量 
级 。 类 似 于 射电 脉冲 星 的 单个 脉冲 但 具有 超过 同一 视线 方向 上 整个 银河 系 星 际 介质 的 预期 总 色散 量 ， 可 作 
为 研究 离子 化 星系 间 介 质 和 自由 电子 空间 分 布 的 宇宙 学 尺度 探 针 。 对 截止 2018 年 6 月 Parkes 望远镜 、 
GBT 望远镜 、Arecibo 电 望 远 镜 、UTMOST 望远镜 和 ASKAP 望远镜 已 探测 到 的 52 例 FRB 事件 进行 了 观 
测量 统计 分 析 。 排 除 来 自 银河 系 星际 介质 贡献 的 色散 量 后 ，FRB 平均 色散 量 约 为 584.521pccm-3 ,这 暗示 
FRB 来 自 河 外 。 通 过 最 佳 拟 合 估算 红 移 窜 律 分 布 谱 dN/dFobs=4.14+1.30xF7L12055sky-1 day-1。FAST 望 
远 镜 使 用 19 波束 接收 机 探测 灵敏 度 后 ,FAST 多 科学 目标 同时 扫描 巡天 ( CRAFTS ) 中 估计 经 过 一 年 时 间 
的 科学 运行 , FAST 能 发 现 大 约 10 个 FRB ,将 有 效 的 扩大 样本 ， 为 FRB 研究 提供 重要 信息 。 
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Fast Radio Bursts Statistics 


Abstract:Fast Radio Burst was first discovered in pulsar search data in 2007, it is intense radio 
flashes, can be a powerful probe for studying the ionized intergalactic medium and the spatial 
distribution of free electrons. After subtracting the contribution from the interstellar medium of our 
Galaxy, Their mean dispersion measure is found to be —584.521pc cm^, supporting their being 
from cosmological origin. The  law-distribution can be described  asdN/dFobs-4.14 + 
1.30xFobs sky! day'! Based on this intensity distribution function, CRAFTS will be able to detect about 10 
FRBs for every a year of observation time. 
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1、 引 言 


2007 年 , 一 类 新 的 射电 暴 首 次 在 脉冲 星 搜索 数据 中 被 发 现 "。 这 类 射电 暴 类 
似 于 射电 脉冲 星 ni 视线 方向 上 整个 银河 系 星际 介质 的 预 
期 总 色散 (Lorimer et al，2007)。 因 只 有 一 个 事例 ， 当 时 并 没有 引起 太 多 的 重 
视 。 直 到 2013 ee em | (Thornton et al. 2013)， 这 类 
爆发 被 命名 为 快速 射电 爆 (Fast Radio Burst, FRB), Thornton 等 人 (2013) 
首次 在 《science》 上 确认 一 个 快速 射电 上 暴 位 置 、 窒 主星 系 以 及 红 移 的 研究 结果 
7, 该 快速 射电 暴 精确 位 置 和 宿主 星系 的 确定 ， 是 继 2016 年 初 LIG0 引力 波 探 测 
之 后 天 文 和 物理 界 最 重要 的 发 现 , 重要 性 在 于 确认 了 射电 观测 现象 和 宇宙 学 尺度 
天 体 物理 起 源 的 直接 连接 。 主 流 学 界 认为 其 来 自 于 宇宙 学 距离 ”， 相 关 研 究 成 为 
射电 天 文 领域 的 热点 。 提 出 的 模型 包括 伽 马 暴 、 黑 洞 并 合 、 巨 脉冲 以 及 地 球 范 围 
内 的 电磁 效应 等 "。 

到 目前 为 共 发 表 了 约 52 个 FRB 的 观测 数据 (FRBCAT, Petroff et al. 2018). 

其 中 22 个 是 Parkes 望远镜 在 1400MHz 波段 发 现 的 ,2 个 是 Arecibo 望远镜 在 1400 
MHz 发 现 的 ，1 是 GBT 望远镜 在 800 MHz 发 现 的 ，4 个 是 UTMOST 望远镜 在 843MHz 
发 现 的 ,其 余 是 ASKAP 望远镜 在 1300MHz 发 现 的 。 扣 除 银河 系 星际 介质 的 贡献 后 ， 
52 个 FRB 平均 的 色散 量 为 584. 521pe cm”， 这 暗示 FRB 来 自 河 外 。 这 些 FRB 

E 量 输出 跨越 了 三 个 量 级 ， 平 均值 约 为 10 ergs， 而 FRB 的 持续 时 间 通 常 为 毫 
秒 量 级 ， 通 常 观测 时 间 持 续 在 I1I-10ms" (Katz et al. 2016)。 由 于 FRB 输出 能 量 
高 ， 持 续 时 间 短 ， 所 以 光度 巨大 。 另 一 方面 持续 时 间 短 也 说 明 这 些 源 非 常 致密 。 
iL 我 们 获得 的 FRB 的 信息 有 限 ， 进 一 步 的 观测 限制 ”是 这 一 领域 的 突破 
方向 (Emily Petroffet al 2017). 

快速 射电 暴 信号 的 色散 延迟 超过 在 视线 方向 上 银河 系 的 色散 量 总 和 。 其 宇宙 
学 尺度 的 河 外 起 源 和 高 射电 流量 密度 的 观测 特征 , 使 快速 射电 暴 的 到 达 时 间 (TOA) 
可 成 为 基本 物理 量 精密 测量 的 灵敏 探 针 。 比 如 ”” (Wei et al. 2015 PRL, 115, 
261101; Keane et al. 2016 Nature, 530, 453): 测定 宇宙 距离 、 测 量 宇宙 中 
的 电子 密度 、 大 尺度 宇宙 学 的 精确 测量 对 暗物质 证 据 、FRB 对 光子 静 质 量 、FRB 源 
属性 、 脉 冲 星 单 脉 冲 对 脉冲 星 辐 射 机 制 、 星 际 内 炬 效应 、 重 子 起 源 和 标准 模型 等 ， 
甚至 检验 爱 因 斯 坦 等 效 原理 乃至 量子 引力 理论 的 强 有 力 工 具 ” 。 随 着 观测 积累 
的 快速 射电 暴 样本 数目 增加 , 甚至 可 通过 样本 色散 延迟 空间 的 多 极 矩 来 分 析 宇 害 
物质 分 布 的 成 团 性 质 ， 测 量 宇宙 大 尺度 结构 。 
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本 文 : 将 在 第 2 节 介 绍 FRB 的 搜寻 过 程 ， 综 述 干扰 消除 、 消 色散 和 特征 信号 
搜索 等 步 又。 第 3 节 对 FRB 样本 进行 了 数据 统计 分 析 和 观测 量 数据 统计 分 析 研 究 
FRB 特性 ， 估 算 FRB 的 强度 分 布 函数 。 第 4 节 讨 论 FAST 对 FRB 的 观测 前 景 。 第 5 
节 进 行 总 结 和 展望 。 


2、FRB 搜寻 的 数据 处 理 


搜寻 核心 思想 是 : 在 有 干扰 信号 的 脉冲 星 了 巡天 或 FRB 巡天 数据 中 , 寻找 具有 
色散 效应 的 短 时 标 、 强 脉冲 性 信号 。 

FRB 信号 数据 量 大 , 内 存 大 和 计算 速度 快 的 epu 或 gpu 保证 FRB 数据 正确 的 
导入 。FRB 数据 记录 的 辐射 流量 密度 Ff， 是 关于 时 间 t、 频 率 f 和 偏振 p 的 多 元 
函数 : 


F-F (t, f, p) (1) 

FRB 信号 射频 干扰 消除 处 理 是 后 期 搜寻 处 理 必 备 前 提 。 根据 FAST 的 19 波束 
设备 馈 源 参数 设置 ,波束 间 互 相 会 有 少量 交 义 重合 区 域 。 因 而 ， 望远镜 不 同 波束 
会 顺 次 经 过 同一 射电 源 , 但 对 同一 源 的 观测 存在 时 间 差 。 为 此 ， 可 利用 多 波束 间 
数据 和 应 用 数据 库 技术 开展 交叉 检验 ,识别 射电 源 并 改善 源 定位 精度 ， 准 确 提 取 
RFI 信息 。 有 具体 算法 采用 多 波束 交叉 相关 和 梯度 下 降 检测 技术 结合 算法 : 多 波束 
交叉 相关 算法 对 两 个 同一 观测 时 间 点 的 不 同 波束 进行 时 间 位 移 的 互相 关 , 得 出 两 
个 矩阵 在 时 间 差 上 的 互相 关系 数 , 并 进行 波束 间 遍 历 和 标准 化 。 每 个 波束 接收 的 
昌 射 功率 谱 可 由 源 与 波束 距离 加 权 平 均 计算 , 再 使 用 梯度 下 降 最 优化 算法 搜索 源 
的 可 能 发 生 位 置 , 最 后 通过 反 卷 积 波 束 形 状 解构 出 多 波束 履 盖 天 区 的 源 强度 分 布 
或 侦 测 获 取 RFI 所 处 频段 分 布 信息 。 

对 任 一 路 偏振 ， 巡 天 数据 以 一 个 二 维 数组 的 形式 存放 于 psrfits 或 者 
fits(The Flexible Image Transport System (FITS) 文 件 中 。 在 基本 的 搜寻 过 
程 中 ， 所 需 完成 的 工作 就 是 在 T( 时 间 )- DM( 色 散 量 ) 组 成 的 两 维 空间 中 搜索 单 脉 
冲 强 流 量 信号 。 当 脉冲 星 信 号 传播 过 程 中 , 会 受到 等 离子 体 中 的 自由 电子 相互 作 
用 导致 群 速度 由 不 同 频率 决定 到 达 早晚 ,低频 信和 号 出 现 延 迟 ， 而 高 频 信 和 号 的 脉冲 
辐射 较 早 抵达 。 故 其 信号 的 到 达 时 间 延 迟 和 频率 之 间 成 平方 反比 关系 。 

Atccv ^ (2) 
不 同 频率 的 电磁 波 在 介质 中 的 传播 速度 不 同 ， 对 于 不 同 频率 f. 信和 号 到 达 延 迟 时 
间 t 为 


ES 
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t-DXTT (3) 


2 


其 中 ，D 为 色散 常数 : D = 


e 
2-Tcme 


= 4.148808 x 10? MHz?pc cms, 色散 量 


DM= 人 medl 表 示 快 速射 电 暴 的 距离 ， 从 D R MARAR 


5 DM = 1.205 x 10 "cm pc.tsonp f^ / ^ f 
消 色 散 完成 后 得 到 相应 的 数据 , 接着 针对 每 1 个 DM 值 产生 的 数据 进行 判别 。 
在 检测 器 中 , 预先 排除 已 知 的 脉冲 星 干扰 并 设置 合理 的 信 品 比 和 时 间 阔 值 参 数 ， 


可 以 得 到 疑似 的 候选 值 ， 再 对 候选 体 筛 选 找到 FRB. 


3、 样 本 和 统计 分 析 


HAEA 


(4) 


综合 当前 澳大利亚 的 帕克 斯 (Parkes ) 望远镜 、 美 国 的 绿 岸 (GBT) 望 远 
镜 、 美 国 的 阿 雷 西 博 电 (arecibo) Hurt, UTMOST 望远镜 (M. Caleet al 2016) 
和 ASKAP 望远镜 共有 52 个 快速 射电 暴 。 其 中 有 一 个 是 重复 的 快速 射电 暴 
FRB121102, 这 个 快速 射电 暴 在 2012 年 被 美国 阿 雷 西 博 射电 望远镜 第 一 次 探测 


到 ， 接 着 在 2015 年 和 2016 年 再 次 被 GBT 和 美国 的 甚大 天 线 聊 
到 爆发 , Tendulkar et al "第 一 次 精确 地 测 出 该 快速 射电 暴 的 红 移 和 宿主 星 


EMIF 25 次 探测 


系 ， 确 认 了 快速 射电 暴 的 宇宙 学 起 源 。 其 余 都 是 非 重 复 的 快 射 射电 暴 ，52 个 
快速 射电 暴 参 数 如 表 1 所 示 。52 个 FRB 最 小 DM 值 为 114. 1, 最 大 的 高 达 2596. 1; 


其 中 Fobs=SpeakXW，DMExcess=DM-DMGalaxy: 


Table 1.1. 截 止 到 2018 年 6 月 己 观测 到 的 52 个 快速 射电 暴 样 木 
FRB telescope | b DM DMGalaxy DMExcess SNR  Wobs Speak Fobs 21 z2 23 
(degree) (pc:cm^3) (pc:cm^3) (pc:cm^3) (ms) (Jy) (Jy:ms) 
FRB010125 parkes 356.641 -20.0206 790+3 110 680+3 17 9.4 0.3 2.82 1.001 0.931 0.861 
FRB010621 parkes 25.434 -4.00381 745+10 523 222+10 16.3 7 0.41 2.87 0.594 0.524 0.454 
FRB010724 parkes 300.653 -41.8051 375 44.58 330.42 23 5 30 150 0.394 0.324 0.254 
FRB090625 parkes 226.444 -60.0303 899.55+0.01 31.69 867.86+0.01 30 1.92 1.14 2.1888 1.226 4.155 1.085 
FRB110214 parkes 290.7 -66.6 168.9+0.5 31.1 137.8x0.5 13 1.9 27 51.3 0.207 0.137 0.067 
FRB110220 parkes 50.829 -54.7663 944.38+0.05 34.77 909.61+0.05 49 5.6 1.3 7.28 1.29 1.22 1.15 
FRB110523 GBT 56.12 -37.82 623.3*0.06 43.52 579.78+0.06 42 1.73 0.6 1.038 0.828 0.758 0.687 
FRB110626 parkes 355.862 -41.7522 7230.3 47.46 675.54+0.3 11 1.4 0.4 0.56 0.967 0.897 0.826 
FRB110703 parkes 80.9978 -59.0191 1103.6x0.7 32.33 1071.27x0.7 16 4.3 0.5 2.15 1.515 1.45 1.375 
FRB120127 parkes 49.2871 -66.2037 553.3*0.3 31.82 521.48x0.3 11 1.1 0.5 0.55 0.747 0.677 0.607 
FRB121002 parkes 308.22 -26.2647 1629.18x0.02 74.27 1554.91 16 5.44 0.43 2.3392 2.199 2.129 2.059 
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0.02 


FRB121102 arecibo 174.95 -0.225138 557+2 188 369+2 14 3 0.4 1.2 0.32 0.32 0.32 

FRB130626 parkes 7.45003 27.4203 952.4£0.1 66.87 885.53+0.1 21 1.98 0.74 1.4652 1.244 1.174 1.104 

FRB130628 parkes 225.955 30.6556 469.88+0.01 52.58 417.3+0.01 29 0.64 1.91 1.2224 0.593 0.523 0.453 

FRB130729 parkes 324.788 54.7446 86122 31 83042 14 15.61 0.22 3.4342 1.172 1.101 1.031 

FRB131104 parkes 260.55 -21.9253 T7941 71.1 707.9+1 30 2.08 1.12 2.3296 0.783 0.713 0.643 

FRB140514 parkes 50.8413 -54.612 562.7£0.6 34.9 527.8+0.6 16 2.8 0.471 1.3188 0.755 0.685 0.615 

FRB141113 arecibo 191.9 0.36 400.3 188 212.3 84 2 0.039 0.078 0.145 0.075 0.005 

FRB150215 parkes 24.6628 5.28092 1105.6x0.8 427.2 678.4+0.8 19 2.88 0.7 2.016 1.135 1.065 0.994 

FRB150418 parkes 232.665 -3.2348 776.2+0.5 188.5 587.7+0.5 39 0.8 22 1.76 0.632 0.562 0.492 

FRB150610 parkes 278 16.5 1593.9x0.6 122 1471.9x0.6 18 2 0.7 1.4 2.062 1.992 1.922 

FRB150807 parkes 333.892 -53.5959 266.5+0.1 36.9 229.6+0.1 0 0.35 128 44.8 0.338 | 0.268 0.198 

FRB151206 parkes 32.6 -8.5 1909.8x0.6 160 1749.8x0.6 10 3 0.3 0.9 2.4452 2.382 2.311 

FRB151230 parkes 239 34.8 960.4+0.5 38 922.4+0.5 17 44 0.42 1.848 1.204 1.223 1.153 

FRB160102 parkes 18.9 -60.8 2596.1+0.3 13 2583.1+0.3 16 3.4 0.5 1:7 3.609 3.539 3.469 

FRB160317 UTMOST 246.05 -0.99 1165+11 319.6 845.4+11 13 21 3 63 1.08 1.01 0.94 

e ) FRB160410 UTMOST 220.36 27.19 278+3 57.7 220.3+3 13 4 7 28 0.31 0.24 0.17 
£3 FRB160608 UTMOST 254.11 -9.54 682+7 238.3 443.747 12 9 4.3 38.7 0.521 0.451 0.381 
> FRB170107 ASKAP 266 54.1 609.5+0.5 35 574.5+0.5 16 2.4 24.1 57.84 0.819 0.749 0.679 
" FRB170416 ASKAP 337.6 -50 523.210.2 40 483.2+0.2 13 5 19.4 97 0.695 0.625 0.555 
FRB170428 ASKAP 359.2 -49.9 991.7+0.9 40 951.7+0.9 10.5 4.4 T.T 33.88 1.852 1.282 1.212 

FRB170707 ASKAP 269.1 -50.5 235.2+0.6 36 199.2+0.6 9.5 3.5 14.8 51.8 0.292 0.222 0.152 

. - FRB170712 ASKAP 329.3 -51.6 312.79+0.07 38 274.79+0.07 12.7 1.4 37.8 52.92 0.401 0.331 0.261 
- FRB170827 UTMOST 303.2 -51.7 176.8+0.04 37 139.8+0.04 110 0.4 60 24 0.211 0.141 0.071 
FRB170906 ASKAP 34.2 -49.5 390.3+0.4 39 351.3+0.4 17 2.5 29.6 74 0.51 0.44 0.37 

FRB171003 ASKAP 283.4 46.3 463.2+1.2 40 423.2+1.2 13.8 2 40.5 81 0.6 0.53 0.46 

FRB171004 ASKAP 282.2 48.9 304+0.3 38 266+0.3 10.9 2 22 44 0.38 0.31 0.24 

i FRB171019 ASKAP 52.5 -49.3 460.8+1.1 37 423.8+1.1 23.4 5.4 40.5 218.7 0.609 0.539 0.469 
: FRB171020 ASKAP 29.3 -51.3 114.1+0.2 38 76.1+0.2 19.5 3.2 117.6 376.32 0.124 0.054 0.01 
FRB171116 ASKAP 205 -49.8 618.5+0.5 36 582.5+0.5 11.8 3.2 19.6 62.72 0.815 0.745 0.674 

FRB171213 ASKAP 200.6 -48.3 158.6+0.2 36 122.6+0.2 25.1 1.5 88.6 132.9 0.167 0.097 0.027 

FRB171216 ASKAP 273.9 -48.4 203.1+0.5 37 166.1+0.5 8 1.9 21 39.9 0.245 0.174 0.104 

FRB180110 ASKAP 7.8 -51.9 715.7402 38 677.7+0.2 35.6 3.2 128.1 409.92 0.967 0.897 0.827 

FRB180119 ASKAP 199.5 -50.4 402.7x0.7 36 366.7x0.7 15.9 2.7 40.7 109.89 0.512 0441 0.371 

FRB180128.0 ASKAP 326.7 522 441.4x0.2 32 409.4x0.2 12.4 2.9 17.5 50.75 0.582 0.512 0.441 

FRB180128.2 ASKAP 327.8 -48.6 495.9x0.7 40 455.9x0.7 9.6 2.3 28.7 66.01 0.656 | 0.586 0.515 

FRB180130 ASKAP 5.9 -51.8 343.5+0.4 39 304.5+0.4 10.2 4.1 23.1 94.71 0.445 0.375 0.305 

FRB180131 ASKAP 0.6 -50.7 657.7+0.5 40 617.7+0.5 13.8 4.5 222 99.9 0.884 0.814 0.744 

FRB180212 ASKAP 338.3 50 167.5x0.5 33 134.5x0.5 18.3 1.81 53 95.93 0.196 0.126 0.056 

FRB180315 ASKAP 13.2 -20.9 479x0.4 36 443+0.4 10.4 2.4 23.3 55.92 0.572 0.502 0.432 

FRB180324 ASKAP 245.2 -20.5 431x0.4 70 3610.4 9.8 4.3 16.5 70.95 0.504 0.434 0.364 


基金 项 目 : 国家 自然 科学 基金 ( U1631132 ) 
收 稿 日 期 : 2018 - 11 - 21; 
作者 简介 : RUR, 女 ,硕士 .研究 方向 : FAST 数据 处 理 研究 Email: gznuyqy@163.com 


FRB180525 ASKAP 349 50.7 388.1+0.3 31 357.1+0.3 27.4 3.8 78.9 299.82 0.506 


0.436 


0.365 


表 1 截止 2018 年 3 月 观测 到 的 52 个 快速 射电 暴 样本 "第 1 列 提供 了 FRB 名称， 第 2 列 指明 观测 的 所 属 望 

远 镜 。 第 3 列 和 第 4 列 分 布 给 出 各 个 快速 射电 暴 的 银 经 C 和 银 纬 (b); 第 5,6 和 7 列 分 别 显示 了 观察 

到 的 快速 射电 暴 的 色散 量 (DM)， 来 自 银河 系 所 贡献 的 色散 量 (DMeus*) MAMRE (Mees); 第 8 列 

ERRE: 第 9 列 给 出 了 观察 到 的 观测 脉冲 宽度 (Wa. 28 10 列 是 FRB 的 峰值 流量 密度 (Sea); 第 11 RE 

以 Jyms 为 单位 列 出 观测 流量 Fo); 第 12, 13 和 14 列 是 分 别 在 Di。 为 0、50 和 100 时 的 红 移 值 。 注 红 移 

计算 通过 YMW16 包含 了 一 个 星系 间 介 质 中 自由 电子 密度 的 模型 ，FRB12110 不 存在 该 假设 ， 其 红 移 已 明确 。 
ik: 表格 源 自 网 站 astronomy. swin. edu. au/pulsar/frbcat 

Tablel Column 1 shows the FRB names, Column 2 indicates the telescope, The longitude and 
latitude of each FRB are given in Column 3, 4. Column 5, 6, and 7 provide the observed DMs of FRBs, 
the DM contributions from the Galaxy and the DMExcess, respectively. SNR is presented in Column 
8. Column 9 provides the observed width duration of the corresponding ratio pulse (Wobs). Column 
10 is the observed peak flux density(Speak) of each FRB.Column 11 tabulates the observed 
fluences (Fobs). Column 5, 6, and 7 provideredshift values at DMhost 0, 50, and 100, respectively. 

Data from http://y. swin. edu. au/pulsar/frbcat/ (Petroff et al. 2018). 

KERREN ERAT HR ZB) a EA 2. 6063* $) 83.2938? , FRBs 在 各 个 天 
"31898 REZ118] E^ BUSERHSE. MZAA AEE Keith W 2014). 。 其 中 ， 存 
在 有 FRBs 空间 分 布 位 置 十 分 接近 如 FRBI40514 和 FRB110220， 两 者 之 间 位 置 仅 
相差 9 角 分 左右 ， 两 者 DM 相差 380pc cm ”。 对 此 Piro&Burke-Spolo 认为 其 是 同 
一 个 FRB, 起 源 于 超新星 遗迹 中 的 中 子 星 ”， 由 于 超新星 喷 流 膨胀 会 导致 天 体 源 
附近 环境 发 生变 化 ， 导 致 色散 量 发 生 改变 ”“。 此 外 比较 相同 纬度 FRB 和 脉冲 星 的 
色散 量 时 发 现 FRB 的 色散 量 显著 偏 高 ,超出 银河 系 所 贡献 的 色散 量 的 10 一 20 倍 ， 
该 情况 被 认为 主要 是 星际 介质 中 自由 电子 所 导致 。 


3. 1 数据 分 析 
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图 1DM 和 DMExcess 的 统计 直方 图 


Fig.1statistical histogram of DM and DMExcess 
对 FRB 色散 量 进一步 观察 ，DM 的 函数 峰值 为 664. 5246pccm ， 而 DMExcess 
峰值 为 584. b21pccm ^. Kr DMExcess 占 DM 的 87. 9596, 支持 了 FRB 是 天 文学 起 
源 。 如 图 1 里 图 a 和 图 b 所 示 的 DM 和 DMExcess 的 直方 图 的 拟 合 曲线 呈 瑞 利 分 布 ， 
HETA 1 的 瑞 利 分 布 情况 。 根 据 瑞 利 公式 


—e 257 
y=} ne ,x>0 (5) 
0 ,xx0 
推导 出 DM 与 FRB 数量 的 关系 ， 公 式 如 下 : 
EE 
N=DMe 2 (6) 
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图 2 分 别 给 出 了 观测 流量 密度 、 峰 值 流量 密度 与 DMexeess 的 相关 性 ， 峰 值 流 量 密 度 、 对 应 
色散 量 与 脉冲 宽度 的 相关 性 


Fig2The correlation between Fobs, Speak and DMExcess respectively. The 
correlation between Speak DMExcess and Wobs respectively 

注 : 误差 棒 的 取 值 源 自 文献 

图 上 Ca). (b) 显示 Fobs 和 Speak 与 DMes 没 有 任何 明显 的 相关 性 ， 

出 乎 理论 上 的 结果 ， 可 能 存在 两 方面 原因 : 1. 星系 际 IM 介质 的 电子 密度 
分 布 不 清 ， 模 型 依赖 严重 ， 需 要 更 多 样本 和 多 波段 观测 约束 和 提高 IGM 电子 密 
度 模 型 的 真实 性 ，2. FRB 能 量 输出 跨越 数 个 量 级 ， 不 能 完全 体现 1/r 的 距离 相 
关 性 也 是 正常 的 。 根 据 图 内 容 大 致 得 到 显示 Fobs 和 Speak 与 DMsweess 呈现 负 相 
KRR. (c) 说 明 部 分 较 强 的 流量 密度 的 快速 射电 暴 有 具有 较 罕 的 脉冲 宽度 ， 
这 些 符合 脉冲 星 的 巨型 脉冲 的 相关 性 质 ，Geng&Huang 提出 ， 小 行星 与 中 子 星 
的 碰撞 接触 面积 小 时 , 发 生 的 碰撞 更 激烈 , 导致 流量 密度 较 大 快速 射电 其 更 亮 ， 
即 碰 接 接 触 面积 大 时 , 伞 撞 过 程 出 现 延 长 ， 造 成 流量 密度 较 小 快速 射电 其 更 瞳 
“， 故 快速 射电 暴 可 能 起 源 于 小 行星 与 中 子 星 碰撞 。 图 (d) 说 明 脉 冲 宽度 和 
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DMsxeess 具有 正 相 关 关 系 ， 脉 冲 宽度 越 大 ，DMss 越 大 。 
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图 3 图 a 观 测 流量 密度 与 红 移 的 相关 性 ， 图 b 对 应 色散 量 与 红 移 相关 性 
Fig3 panels (a) correlation between Fobs and z. panels(b) correlation between 


DMExcess and z 


采用 了 YMW16 包含 了 一 个 星系 间 介 质 中 自由 电子 密度 的 模型 , 计算 FRB 红 
移 。 无 法 确定 FRB 是 否 源 于 宿主 星系 ， 所 以 DMhost 的 不 同 取 值 作为 红 移 的 统 
计 误 差 ， 选择 的 IGM， 目 的 是 为 避免 脉冲 星 的 DMs 大 于 该 方向 的 最 大 YMW16 模 
型 预测 值 ( 显 然 不 正确 ) 。 图 Ca) 明显 看 到 红 移 和 观测 流量 没有 相关 的 关系 ， 
图 Cb) 红 移 和 DM 成 线性 关系 ， 随 着 DM 的 增 大 红 移 随 着 加 大 ， 通 过 红 移 和 DM 
都 可 大 致 确定 射电 暴 所 在 。 此 外 ， 重 复 射 电 暴 FRB121102 的 重复 爆发 被 确定 了 
天 区 明确 宿主 星系 和 红 移 最 有 可 能 的 值 是 0. 32 (Cordes, Lazio) "^. 研究 学 
者 关于 该 重复 射电 暴 认为 有 可 能 的 前 身 星 物 理 模 型 有 软 伽 马 重复 暴 "、 磁化 中 
子 星 穿 过 小 行星 带 史 以 及 中 子 星 -白矮星 双星 系统 等 ” 。 


3. 2 强度 分 布 


曾 统计 FRB 具有 非 泊 松 重复 ”， 且 现实 爆发 更 加 复杂 ， 从 观察 推断 出 的 事件 
率 取决 于 观察 策略 (Liam Connor 2016)。 当 前 数据 采取 窜 律 函数 的 统计 分 布 情 
况 , 将 一 条 斜率 为 桥 指 数 的 负数 的 直线 , 也 就 是 霸 律 函数 去 是 判定 FRBs 强度 分 布 ， 
助 于 今后 发 现 更 多 的 FRB。 采 取 Tavanid ”对 伽 马 射 电 暴 的 峰值 流量 密度 的 累计 
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分 布 结果 , 因为 在 峰值 流量 密度 较 小 的 一 段 改 分 布 明显 的 偏离 欧 几 里 得 平坦 空间 
的 Nee S. ^ 关系， 然后 改 结论 被 伽 马 暴 红 移 被 精确 测定 后 证 实 是 正确 ， 考 虑 
对 快速 射电 暴 进行 类 似 研究 ， 此 外 Bera et al", Wang&Yu"", Lilongbiao "5&5 
许多 研究 者 将 Fobs 作为 统计 分 析 和 宇宙 学 应 用 的 选择 。 故 将 每 天 在 特定 的 效应 
Fobs 下 出 现 的 FRB 的 实际 数 密 度 去 做 暴 律 分 布 分 析 ， 公 式 如 下 


二 = AF3% CD 


其 中 A 是 以 事件 sky 天 为 单位 的 常 系数 和 a ERRIETA AFI a 都 需要 通 
过 观 紧 来 确定 。 


(b) 


Power-Law Index 


10? 10i 10? 2 4 6 8 10 
Fobs(Jy ms) AF(y ms) 


图 4 Ca) 表示 观测 流量 密度 组 距 为 2. 2 时 的 强度 分 布 函数 , X 轴 的 误差 棒 代 表 组 距 的 大 小 ， 
Y 轴 表 示 每 个 区 间 对 应 的 FRB 数目 ， 其 误差 棒 为 统计 误差 ， 即 区 闻 内 样本 数目 的 平方 根 值 。 
(b) 将 FRBs 按照 不 同 的 Fobs 组 距 划分 后 拟 合 得 到 的 宕 律 指数 a。 

Fig4 panels(a) provides the intensity distribution function when AF-2.2,The error bar on the 
X-axis represents the size of the spacing, and the Y-axis represents the number of FRBS 
corresponding to each interval. The error bar is the statistical error, namely the square root 
value of the number of samples within the interval. (b) the power law exponent a obtained by 


fitting FRBs according to different Fobs group spacing partition. 

因为 当前 FRB 的 总 数量 不 足够 多 ， 所 以 组 距 的 选择 至 关 重 要 。 将 52 个 非 重 
复 快 速射 电 暴 的 峰值 流量 密度 宽度 尝试 各 种 不 同 值 来 研究 其 效果 ， 因 为 FRBs 的 
观测 流量 值 范 围 跨度 较 大 , 所 以 当 组 距 非常 小 CA Fobs<2) 时 , 分 组 的 数目 多 ， 
数据 点 过 于 分 散 ， 拟 合 时 效果 极 差 ， 不 可 用 于 分 析 强 度 分 布 函数 ， 当 组 距 过 大 
(A Fobsz8) 时 ,分 组 误差 较 大 ， 计算 的 结果 没有 意义 ; 中 间 存 在 一 段 组 距 (3 
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XA Fobsx7.3) 时 ， 每 个 分 组 的 数目 跳动 过 大 ， 拟 合 的 线 穿 过 的 数据 太 少 ， 寡 
律 指数 太 大 , 将 其 作为 异常 值 .A Fops 在 2 到 3 以 及 7.3 到 8 之 间 的 范围 内 时 ， 
最 适合 用 来 计算 宫 律 分 布 函 数 的 参数 。 统计 合适 的 组 距 里 的 斜率 , 采用 和 斜率 的 平 


均值 作为 暴 律 指数 a， 根据 图 4 知 a 取 1.12055。 


FLimit R CF) Reference Derived coefficier 
(Jyms) (sky day ) (10°sky day ) 
3.0 10 Thornton et al 2013 5. 61 士 2.04 
0. 35 3. 1X10' Spitler et al 7. 88 +6. 92 
2.0 2.5X10 Keane et al 2015 1.17 土 0. 51 
1.8 1.2X10” Law et al 2015 5.37 土 2. 47 
4.0 4. 4X 103 Rane et al 2016 2. 863-0. 90 
0. 13-5. 9 6. 0X 103 Champion et al 2013 1. 943-1. 27 
0. 03 3.3X10' Long-Biao Liet al 2016 4. 143: 1. 30 
表 2 文献 中 FRBs 的 A 值 


Table 2A values of FRBs in literature 


选择 FAST 设备 来 研究 FRB 的 探测 情况 ， 采 取 A=4. 14 土 1. 30X10'sky day” 
(Long-Biao Liet al 2016)。 故 FRB 的 强度 分 布 函 数 为 : 


二 414 土 1.30 x 103 x Egi sky-1day-! (8) 


4、FAST 对 FRB 的 观测 前 景 

射 天 天 文学 的 发 展 必 备 高 灵敏 度 的 大 口径 射电 望远镜 , 目前 在 大 型 望远镜 巡 
天 设备 有 美国 305m 射电 望远镜 Arecibo 、 德 国 Effelsberg100m 射电 望远镜 、 
澳大利亚 64m 射电 望远镜 Parkes 和 中 国 500m 射电 望远镜 等 ， 其 中 中 国 的 Fast 
是 建 于 贵州 平 塘 的 500 米 口 径 射 电 球面 望远镜 , 利用 天 然 的 喀斯特 洼 坑 作 为 台 址 ， 
数 千 块 可 调节 的 单元 构成 500 米 球 冠 主动 反射 面 ,作为 世界 上 最 大 、 最 灵敏 的 单 
口径 射电 望远镜 ,将 在 中 性 氧 巡 视 、 脉 冲 星 搜 索 、 国 际 VLBI 网 联 测 及 地 外 生命 
搜寻 等 重要 前 沿 领域 取得 突破 9。FAST 天 空 覆盖 的 顶峰 角度 为 26.4” (Li di et 
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al 2018) ,观测 频率 为 70 MHz ~3 GHz (未 来 可 升级 到 8GHz)， 接 收 机 波束 数 为 19 (未 来 焦 
平面 阵列 大 于 100)， 灵 敏 度 (L 波段 ) 为 18K.Jy- ， 系 统 温 度 为 Tsys<*50K， 天 空 覆 盖 范 围 为 
赤 纬 -14" 12' «óx65^ 48' 。 和 Arecib 望远镜 相 比 ，Fast 的 原始 灵敏 度 高 2.5 倍 
Cn, WREE 10 倍 ， 综 合 灵敏 度 高 10 倍 ， 履 盖 的 天 区 是 2 到 3 4%. FAST 的 
接收 器 将 覆盖 低频 (70 到 500MHz) 和 中 频 (0.5 到 3GHz) 频段 。Fast 有 助 于 
发 现 更 多 FRB， 对 于 更 深 一 步 的 FRB 研究 也 会 大 有 帮助 。 在 探测 毫秒 脉冲 星 、 肪 
冲 双 星 、 双 脉冲 星 和 银河 系 外 的 脉冲 星 等 方面 更 具 优 势 。FAST 还 具有 发 现 脉 冲 星 黑 洞 体系 
的 科学 潜力 ， 这 将 有 助 于 了 解 黑洞 的 精确 信息 。 

FAST 的 19 波束 的 优势 和 特点 ， 有 助 于 开展 多 科学 目标 同时 扫描 巡天 : 脉 
冲 星 搜索 、 气 体 星系 搜索 、 银 河 系 中 性 氧 成 像 、 快 速射 电源 的 搜索 。 世 界 上 其 他 
射电 望远镜 ,还 没有 这 四 大 目标 同时 进行 ,是 世界 上 没有 前 例 的 开创 。 构 建 出 比 
较 有 效 的 扫描 巡天 模式 ， 提 供给 全 国 天 文学 者 ， 产 生 系统 性 的 发 现 和 重大 突破 。 
多 科学 目标 同时 扫描 巡天 测量 FRB 还 具备 实时 触发 机 制 ， 可 将 基带 数据 再 1 分 
钟 内 完成 保存 在 环形 缓冲 区 中 。 具 有 前 所 未 有 的 快速 、 最 小 的 系统 控制 复杂 性 和 
RFI 降低 保证 获得 最 佳 系统 性 能 

通过 预期 PAST 的 观测 限定 FRB 的 发 生 条 件 ， 给 出 FRB 发 生 概率 的 估计 ， 在 
不 考虑 台 址 射电 环境 干扰 的 理想 条 件 下 , 根据 FRB 事件 的 全 天 爆发 率 ， 和 FAST 
采用 19 波束 接收 机 探测 灵敏 度 、FAST 可 见 天 区 范围 C 41000 平方 度 ) Aa 
尺寸 (19 XRM 0.15 PYJE), FAST 的 灵敏 度 或 流量 密度 下 限 可 通过 下 式 
估算 (zhang et al. 2015) 2， 


Suo HA PCS EC EE 


T sys 
其 中 , T ERA, ATER, SNR 为 信 噪 比 ,通常 取 10 作为 可 靠 
的 FRB 检测 (Champion et al. 2016), Az 是 观测 带宽 ，A. 是 有 效 面积 。 由 (Li et 
al. 2017) 知 ， FAST 的 观测 流 LE 量 限 制 为 Fiinit=Siimt X ATE 0. 03Jy ms. 


R-Af"* Fedre C0) 


Flimit 
其 中 ，Flimit 是 射电 望远镜 的 灵敏 度 ，Fmax A Dos] re a8 JU 5E BS ET 


4 0.15deg? _1 _ E 
Rast ~3.578 + 1.30 x 10” x 41253deg 7 69 = 0.13 day ` (11) 


根据 上 式 , 可 计算 出 FAST 在 1000 小 时 内 可 探测 到 约 6 个 FRB,，FAST 19 波 
束 每 日 巡天 8 小 时 将 采集 到 约 5TB 量 级 的 数据 。 全 年 可 工作 200 余天 一 年 的 


yu 


(9) 
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观测 时 间 可 以 探测 到 约 个 10 FRB. 


5、 总 绪 与 展望 


本 文 统计 分 析 了 Parkes 望远镜 、GBT 望远镜 、Arecibo 电 望 远 镜 、UTMOST 
望远镜 和 ASKAP 望远镜 探测 到 的 截止 2018 年 的 52 个 FRB 关键 参数 。FRB 空间 分 
布 呈现 各 向 同性 , DM 范围 从 114. 1 到 2596. 1 pccm-3, 峰值 为 664. 5246 pecm-3, 
而 DMExcess 峰值 为 584. 521 pccm-3。 而 占 DM 总 量 的 约 87. 95%。 表 明 Speak, 
Fobs 与 DMExcess 没有 任何 明显 的 相关 性 ， 存 在 有 DM 和 红 移 呈现 正比 关系 。 

将 观测 流量 作为 FRB 强度 的 因素 之 一 , 计算 快速 射电 暴 事件 发 生 率 , 计算 出 
震 律 分 布 函数 的 指数 a 约 为 1. 12055， 强 度 分 布 函 数 dN/dFobs-4.14 + 1.30 x 
103 x Fp sky" day 。 其 中 存在 有 几 个 因素 误差 。 第 1， 目前 观察 到 的 FRB 
总 数 仍然 非常 有 限 。 它 可 能 会 导致 测量 的 圭 律 指数 有 较 大 的 波动 。 至 少 需要 约 
100FRB, 才能 提取 有 关 FRB 物理 性 质 的 结论 性 信息 。 第 2, FRB 不 是 理想 的 标准 。 
第 三 , FRB 可 能 不 是 同 质 的 光源 ,， 共 移动 密度 或 其 亮度 可 能 会 在 太空 中 演变 。FRB 
事件 实际 上 可 能 实际 上 发 生 在 天 空中 , 但 由 于 当前 的 技术 限制 , 我 们 无 法 记录 或 
找到 它们 ， 未 来 将 需要 更 多 新 的 FRB 样品 。 

综 上 所 诉 , 计算 出 的 强度 分 布 函 数 还 只 是 一 个 中 期 阶段 的 结果 。 以 后 需要 更 
多 样本 来 更 准确 地 确定 寡 律 指数 。FAST 望远镜 具有 大 口径 面积 和 多 波束 探测 ， 
对 FRB 的 探测 有 很 大 帮助 ， 有 助 于 增加 发 现 FRB 样本 。 更 重要 的 是 ，FAST 可 以 
在 很 宽 的 频率 范围 内 工作 , 并 且 可 以 为 FRB 提供 详细 的 频谱 信息 , 通过 预期 FAST 
的 观测 限定 FRB 的 发 生 条 件 ， 给 出 FRB 发 生 概 率 的 估计 ， 计 算出 FAST 在 1000 
小 时 内 可 探测 到 约 6 个 FRB， 一 年 的 工作 时 间 内 探测 到 10 个 FR, 有 望 成 为 该 领 
域 的 强大 工具 。 
致谢 
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